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Quantum Rabi model is exactly solved by employing the parameter-dependent unitary transfor- 
mation method in the Bargmann space. The complete energy spectrum consists of two double-fold 
sub-energy spectra. The eigenvalue is determined by the parameter in the unitary transformation, 
which satisfies a highly nonlinear equation. Such the energy spectrum completely coincides with 
that obtained in the occupation number representation [D. Zhang, chinaXiv:201708.00168]. 
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Quantum Rabi model describes the response of a two- 
level atom to an applied bosonic field [1]. The simplest 
light-matter interacting model has been widely applied 
in atomic physics [2], quantum optics [3], quantum infor- 
mation and quantum technology [4], and etc.. Its Hamil- 
tonian reads 


H = wata + gla + a)oz + Aoz + €0n, (1) 


where c, and o, are the Pauli matrices, at and a are the 
creation and annihilation operators for the single boson- 
ic mode with frequency w, respectively, 2A is the level 
splitting for the two-level system, g is the light-matter 
interaction strength, and the last term eo, is the driving 
term, which leads to tunnelling between the two levels. 
The ratio between g and w divides three experimental 
regimes: (a) the weak coupling regime (g/w <~ 0.1), (b) 
the ultrastrong coupling regime (~ 0.1 < g/w <~ 1.0), 
and (c) the deep strong coupling regime (g/w >~ 1.0). 
We note that most experiments have been performed in 
the the weak coupling regime, where the quantum Rabi 
model (1) with e = 0 is equivalent to the exactly solvable 
Jaynes-Cummings model in the rotating-wave approxi- 
mation [5]. Recently more and more attention has been 
focused on the strong regimes due to their fundamental 
characteristics and the potential applications in quantum 
devices [4]. 

It is known that exact solution of quantum Rabi model 
(1) has been a long-standing problem in Physics. Such 
the analytical solution is expected to explore accurately 
the light-matter interaction from weak to extreme strong 
regime. In Ref. [6], Braak presented an analytical solu- 
tion of the Rabi model by employing the representation 
of bosonic operators in the Bargmann space of analytical 
functions. The energy spectrum consists of two parts, 
i.e. the regular and the exceptional spectrum. Howev- 
er, it has been proved that the Braak’s analytical solu- 
tion is incorrect due to the derivation error in solving the 
time-independent Schrodinger equation in the positive 
and negative parity parts [7]. In Ref. [8], we obtained 
the exact solution of the Hamiltonian (1) in the whole 
range of the physical parameters by using the parameter- 


dependent unitary transformation technique in the occu- 
pation number representation. Such a direct and power- 
ful approach has also been used to solve successfully the 
complex two-dimensional electron gas in the presence of 
both Rashba and Dresselhaus spin-orbit interactions un- 
der a perpendicular magnetic field [9,10]. The complete 
energy spectrum of the quantum Rabi model consists of 
two double-fold sub-energy spectra. 

In this work, we reinvestigate the eigenvalue problem 
for the Hamiltonian (1) in the Bargmann space, where 
the bosonic creation and anihilation operators in terms 
of a complex variable z can be transformed as at > z 
and a + 0/0z, respectively. In this representation, the 
state W(z) can be normalized according to 


< |Y >= L [dzdze U* (2) (2) =1. (2) 


We shall see below that the energy spectrum is complete- 
ly consistent with that in the occupation number repre- 
sentation [8]. 

In the Bargmann space, the two-component eigenstate 
of the Hamiltonian (1) for the nth energy level with quan- 
tum number s has the general form 
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where s = +1, Ans is a real parameter in the unitary ma- 
trix to be determined below by requiring the coefficients 
A! and B?* to be nonzero. When i + Hoo, A?* > 0 
and BP* — 0. Substituting the eigenfunction (3) into 
the eigen-equation HW,,, = EnsYns and requiring the 
coefficients of zê to be zero, we obtain the infinite system 
of homogeneous linear equations with the variables Aps 
and B!” 
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where 4 = 0,1,2,:::,00, A25 = BS = 0 form < 0. 


It seems very hard to solve Eqs. (4) and (5). Howev- 
er, we can acquire all the physical eigenvalues and the 
corresponding eigenfunctions by using a special trick. 

The sub-energy spectrum I. In order to get the analyt- 
ical solution of the Hamiltonian (1) in the whole param- 
eter space, we first choose 
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in the eigen-equations (4) and (5). Eqs. (6) and (7) 
come from the vanishing of the two terms about AP$: 
and Br’ in Eq. (4) with i = n+ 1 and Eq.(5) with i = n, 
respectively. Such a choice is based on the observation of 
exact solution of the Hamiltonian (1) for the nth energy 
level with quantum number s when g = 0. Then the non- 
zero eigenfunction associated with the eigenvalue Ens is 
solely fixed by requiring 


[2AAns + e(1 — A2 /)]B!? — 2gAns(n + YAY, = 0, (8) 
or 
2gAns B75 + [2ZAAns + e(1 — AZ NANO. (9) 


We solve the homogenous linear equations (6) and (7) 
about Ay, and B? by vanishing the coefficient deter- 
minant. Then the eigenvalue for the nth eigenstate with 
s has the analytical expression 


Ens = (n + iw + sEns, 


NG-AZ ,)—2eAns 1—4}, 
Ens = y (8 + Te)? + (n + DA AE) 


Note that the quasiparticle energy E,,, must be larger 
than zero. Combining Eqs. (6)-(10), the parameter Ans 
is determined by the constraint 


e(1+ A?) — 2Ans(Ens — wn) = 0, (11) 
or 
e(1+ A?,) +2Ans[Ens —w(n+1)}=0. (12) 


After analysing carefully, we find that Eq. (11) with 
s = —1(1) coincides with Eq. (12) with s = 1(-1). So 
we have 


e(1 + A?) + Ans(20Ens — w) = 0, (13) 


FIG. 1: (Color online) The low-lying energy levels of the sub- 
energy spectrum I in unit of w as a function of the coupling 
parameter g when A = 0.4w and e = 0. The solid lines denote 
n = 0,1,--- ,5 and s = 1 while the dash lines mean n = 
1,2,---,5 and s = —1. 


where g = +1. It is absolutely surprising that such the 
analytical expression (10) of E,,, with the parameter con- 
straint (13) has been obtained in the occupation number 
representation [8]. 

Obviously, Ans is independent of s in Eq. (13). When 
e- 0, then A,,, = 0. Therefore, the eigenvalue E,,, has 
a simple form 


Ene = (nt sors (G++ (9 


in the absence of the driving term. It is easy to conclude 
from Eq. (13) that the analytical solution (10) of Ens is 
physical if and only if it must satisfy the limit A,,, > 0 
when e- 0. In order to compare conveniently with the 
energy spectrum presented by Braak, here we use the 
physical parameters in Ref. [6]. Fig. 1 plots the low- 
lying energy levels as a function of g at A = 0.4w and 
e = 0. We can see that this energy spectrum possesses 
the level crossings between the neighboring eigenstates, 
which is dramatically different from that in Ref. [6]. 
When e # 0, Ans in Eq. (13) with o = 1 has an w- 
dependent solution. The corresponding eigenvalues Ens 
(n = 0,1,2,--- ,co,s = +1) form the sub-energy spec- 
trum I. In Fig. 2, we show the low-lying energy lev- 
els and the corresponding parameter A,,, of the sub- 
energy spectrum I as a function of g at A = 0.7w and 
€ = 0.2w. When g = 0, the sub-energy spectrum I re- 
covers the exact eigenvalues for the interactionless case, 
i.e. Eng = wn + sv A? + € with Ans = (A— VA? + &)/e. 
Obviously, if € > O, Ans -? 0. We note that another 
w-dependent solution A,, in Eq. (13) with o = —1 is 
nothing but the sub-energy spectrum II (see below). 
For the eigenstate associated with the eigenvalue (10), 


from Eq. (7), we have 
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(15) 
where B?* is a constant to be determined by the nor- 
malized condition (2). The coefficients a!” and 875, pro- 


portional to B”*, are uniquely obtained by the recursion 
relations 
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for i =n+1,n+2,---,+00. Here we have defined 
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where J is the 2 x 2 unit matrix. 
The sub-energy spectrum II. Now we use another choice 
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in the eigen-equations (4) and (5). Eqs. (19) and (20) 
originate in the vanishing of the two terms about AzS 
and Bri, in Eq. (4) with i = n and Eq.(5) with i = 
n + 1, respectively. The corresponding eigenfunction is 
uniquely determined by letting 


[2AAns te(1— AZ )JAn’ +2gAns(n+1)BP4, SO, (21) 
or 
—2gAns AS + [2AAns + €(1— A?,)|BPS, =0, (22) 


Solving Eqs. (19) and (20), we have the eigenvalue for 
the nth eigenstate with s 


Eng = (n+ $)wt sOns, 
NG-A2 ,)—2eAns 1—42, 
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Here Ans satisfies the nonlinear equation 


e(1 + A?) +2Ans(Ens — wn) = 0, 
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FIG. 2: (Color online) The low-lying energy levels of the sub- 
energy spectrum I in unit of w as a function of the coupling 
parameter g when à = 0.7w and e = 0.2w, shown in (a). The 
solid lines denote n = 0,1,:::,5 and s = 1 while the dash 
lines mean n = 1,2,--- ,5 and s = —1. The corresponding 
Ans are displayed in (b). 


or 


e(1+ A?,) — 2Ans|Ens — w(n +1) =0, (25) 
which is derived from Eqs. (19)-(22). Similar to Eqs. 
(11) and (12), Eq. (24) with s = —1(1) is also consistent 
with Eq. (25) with s = 1(—1). Therefore, we obtain 


e(1 + A?) + Ans(2T@ns +w) =0, (26) 


where 7 = +1. We note that the analytical expression 
(23) of Ens and the parameter equation (26) are also 
consistent fully with those presented in the occupation 
number representation [8]. 

Ife = 0, then Ans = 0 from Eq. (26). So the eigenvalue 
(23) also has an explicit expression 
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FIG. 3: (Color online) The low-lying energy levels of the sub- 
energy spectrum II in unit of w as a function of the coupling 
parameter g when A = 0.4w and e = 0. The solid lines denote 
n = 0,1,--- ,5 and s = 1 while the dash lines mean n = 
0,1,:::,5 and s = —1. 


We depict the low-lying energy levels as a function of g 
at A = 0.4w and e = 0 in Fig. 3. 

Ife £0, Ans in Eq. (26) with 7 = 1 also has an w- 
dependent solution. The corresponding eigenvalues con- 
stitute the sub-energy spectrum II. Fig. 4 exhibits the 
low-lying energy levels of the sub-energy spectrum II as 
a function of g at A = 0.7w and e = 0.2w and the corre- 
sponding parameter Ans. When g = O, the sub-energy 
spectrum II also recovers the exact eigenvalues for the in- 
teractionless two-level system, i.e. Ens = wn+svVd2 + €? 
with Ang = (A — VA? + @)/e. We note that after tak- 
ing the transformation Ans > —1/A,., the eigenvalue 
(23) with 7 = 1(-1) becomes the eigenvalue (10) with 
o = —1(1). Therefore, both the sub-energy spectrum I 
and II are double degenerate. 

For the nth eigenstate with s in the sub-energy spec- 
trum II, from Eq. (18), we have 


ns _ (1+4? 
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(28) 
where A?* is the normalized constant. The coefficients 
A!” and B!”, proportional to A?*, obey the same recur- 
sion relations (16) and (17) in the sub-energy spectrum 
I. 

In summary, we have exactly solved the quantum Rabi 
model (1) in the Bargmann space. The complete ener- 
gy spectrum is comprised of two double-fold degenerate 
sub-energy spectrum I and II. Such the exact solution 
can help us to deeply understand the light-matter inter- 
action, especially in strong coupling regimes. Because 
the eigenvalue E,,, has the same analytical expression 
in both the Bargmann space and the occupation num- 
ber representation, this guarantees the correctness of the 
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FIG. 4: (Color online) The low-lying energy levels of the sub- 
energy spectrum II in unit of w as a function of the coupling 
parameter g when à = 0.7w and e = 0.2w, shown in (a). The 
solid lines denote n = 0,1,:::,5 and s = 1 while the dash 
lines mean n = 0,1,---,5 and s = —1. The corresponding 
Ans are displayed in (b). 


exact solution for the quantum Rabi model. 
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